Boundary Green's function for spin-incoherent interacting electrons in one dimension 
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The spin-incoherent regime of one-dimensional electrons has recently been explored using the 
Bethe ansatz and a bosonized path integral approach, revealing that the spin incoherence dramati- 
cally influences the correlations of charge excitations. We here introduce a bosonization scheme for 
strongly interacting electrons, allowing us to generalize the description to account for the presence of 
an open boundary. By calculating the single-electron Green's function we find that the charge sector 

O power-law scaling is highly sensitive to the boundary. Our result allows for a detailed description of 
the crossover between boundary and bulk regimes. We predict that scanning tunneling microscopy 
, on a spin-incoherent system will pick up oscillations in the differential tunneling conductance as 

£SJ ■ a function of the applied voltage V at "intermediate" distances x from a real or a dynamically 

generated boundary. The wavelength of the oscillations, irv c /x, probes the speed v c of the charge 
excitations, and therefore the strength of the electron-electron interaction. 



PACS numbers: 71.10.Pm, 71.27. +a, 73.21.-b 



I. INTRODUCTION 



i 

S: 

The spin-incoherent regime of onc-dimcnsional strongly interacting, very low-density electrons has recently attracted 
' O ■ a lot of interest^. For zero temperature the kinetic energy of an electron can be estimated by the Fermi energy 
Ep — (nhn) 2 /8m. For low densities, n <C a^ 1 , this energy is small compared to the Coulomb potential energy e 2 n/e 
(with e being the dielectric constant and = eh 2 /me 2 the effective Bohr radius of the material). In the limit of 
' low densities the system turns into a Wigner crystal, which — in a classical picture — can be viewed as a system of 
electrons placed equidistantly so as to minimize the potential energy. Quantum fluctuations induce an exponentially 

■ small antiferromagnetic spin exchange J > between the electrons. In the case when J is the smallest energy scale, 
J -C T <C Ep, the spin exchange can no longer support collective spin excitations, and the system is driven to a spin 

■ incoherent regime. The physics of the spin-incoherent regime has been addressed using the Bethe ansata 3 -^ and a 
bosonized path integral approach^. Surprisingly, it was found that the spin incoherence dramatically influences the 
correlations of charge excitations, leading to a power-law decay in the charge sector with an interaction-dependent 
nonunitary exponent. 

In this paper we generalize the description in Ref. [5| to account for the presence of an open boundary. For this 
purpose we introduce a bosonization scheme, valid for strongly interacting electrons, and presented in Sec. |TTJ In 
Sec. IIIII we then derive the single-electron Green's function in the presence of an open boundary condition. Scanning 
tunneling microscopy (STM) as a possible experimental probe of the spin-incoherent regime is discussed in Sec. IIVI 
The last section contains a brief summary of our results. 



T3 

a 

Q II. BOSONIZATION IN THE STRONG COUPLING REGIME 

o ' 

We here introduce a bosonization scheme for strongly interacting electrons, for which the applicability of an ordinary 
bosonization^ becomes questionable. In the strong-coupling regime, the electrons are localized at the lattice sites of 
the Wigner crystal, and the usual procedure of linearizing the spectrum around the Fermi points ±kp and then 
■ expanding the electron fields in left and right movers is no longer justified. In what follows we show that one can 
nonetheless perform an "effective" bosonization, valid for low energies and large distances. 

To begin with, it can be shown that a bosonic Hamiltonian correctly describes the low-energy properties of a 
one-dimensional (ID) Wigner crystal!. To justify this statement we turn to the classical picture and describe a ID 
Wigner crystal as a system of electrons vibrating around their equilibrium positions. This is very similar to the 
classical description of phonons as lattice vibrations. For low energies (long-wavelength limit) the vibrations can be 
described by elasticity theory, with the energy of the system expressed by the Hamiltonian 



H = dx 
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' mns 2 (d x u) 2 
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Here u(x) is the displacement of the medium and p(x) is the momentum density. Density fluctuations are given by 
Sn = —nd x u, and s plays the role of the speed of density waves. 

The classical Hamiltonian in Eq. ([1]) can be straightforwardly quantized by imposing a canonical commutation 
relation between u(x) and p(x), 

[u(x),p(x')]=i5(x-x'). (2) 

The resulting quantum Hamiltonian describes the propagation of density fluctuation in the Wigner crystal, and can 
be written in terms of a bosonic field tp c {x) and its conjugate momentum n c (a;), connected to u(x) and p{x) by 

u ( x ) = i^Vcix), P{x) = 7 =-IL c (x). (3) 



It follows that 

■ji r r 1 

(4) 



K c ii 2 c + ^-{d xVc f 



where 



s, = — , (5) 



with v c the speed of the propagating charge density fluctuations, and with vf — nn/2m the Fermi velocity for 
noninteracting spinful electrons. 

We have thus managed to write the Hamiltonian for a ID Wigner crystal in bosonized form although we did not 
know the bosonization formula which connects electron and boson fields. We next derive such a formula, applying 
a kind of "reverse engineering" to the bosonic description of the Wigner crystal above. First we observe that in the 
J — > limit, electrons behave like spinless fermionic particles, i.e. only charge degrees of freedom survive and spin 
degrees of freedom simply define the huge degeneracy of the ground state. We may thus omit the spin index and write 

*(x)=i/>_(x) + il> + (x), (6) 

where tp-( x ) bk+{ x )] is the part of the electron operator which contains negative [positive] momenta. In analogy with 
ordinary bosonization of spinless fermions^ we introduce two chiral boson fields 4>c- and 4>c+, connected to ip c by 
<p c = 4> c - + 4> c +- We can then write an effective bosonization formula, valid for low energies, 

V27ra 

with a a short-distance cutoff, and where kp and A are to be determined. Using that the density operator is given by 

p=&(x)V(x)=n+^d x (p D , (8) 

with the product of electron fields defined by point splitting^, and where n is the average electron density and 
(y // 2/ir)d x ipc is a fluctuation term [cf. Eqs. ((TJ and j3])], we find that 

kp = 7m, A = 8n. (9) 

We have here used that the fast oscillating nonchiral terms ip±ip T contained in $^'5 effectively average out to zero over 
large distances, and hence can be neglected in the long-wavelength limit. Also note that the doubling of the Fermi 
momentum in Eq. ([9]) is in agreement with our description of spinless fermions, since now only a single fermion can 
occupy a given momentum state. As can be easily verified, the case of noninteracting spinless fermions corresponds 
to K c — 1/2. It follows that in this bosonization scheme interactions are absorbed in two steps: Local interactions 
among the spinful electrons are incorporated in a free Hamiltonian for spinless fermions (with K c = 1/2), while 
long-range interactions renormalize the value of K c (away from K c = 1/2). To see this explicitly, note that for a delta 
function interaction among the electrons — corresponding to the low-energy, long-wavelength limit of the infinite-C/ 
ID Hubbard model^ — the description as noninteracting spinless fermions becomes exact (with K c — 1/2). 
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To obtain a more conventional parametrization (where a unit value of an effective "charge parameter" K corresponds 
to the case of noninteracting spinless fermions) we define a bosonic field ip = y/2(p c with conjugate momentum 
II = n c /\/2- Then the bosonization formula in Eq. (J7J) and the Hamiltonian in Eq. ([3]) take the forms 



V2t 



± 



(10) 



and 



H = 7 ' 



#n 2 + — (0 a¥ >) 



(11) 



respectively, with if = 2K C — 2vf/s = Vp/s (in agreement with the doubling of the Fermi momentum kp). 



III. BOUNDARY GREEN'S FUNCTION 



The bosonization procedure presented in the previous section can easily be adapted to the case when a boundary 
is present. Imposing an open boundary condition (OBC) at the end, x = 0, of a semi-infinite system with x > 0, we 
have that 

*(0) =^_(0)+^ + (0) =0. (12) 

Analytically continuing the chiral fermion operators to negative coordinates, 

1>+(x) = -1>_(-x), (13) 

and then following the standard procedure for an OBC2. — using our Eq. (JTTJJ) — we obtain a bosonization formula 
for strongly interacting spin-incoherent electrons with an OBC, 

1p (ar) 1=3 1 c -iie F x c -i\/iir\cosh('d)<j>- (z) -sinh(fl) jL (-x)] 

\Z2na 

Here e 2 * 5 = if, with = cosh(i9)^>_ — sinh($)0_|_ governed by the free chiral boson Hamiltonian 

#=| fdx[dj-(x)} 2 . (15) 

To calculate the zero-temperature single-electron Green's function in the presence of the OBC we follow the path 
integral approach introduced in Ref. Q. The averaging over spin introduces a factor 2~ N ( X ' X ,r ' in the expression for 
the Green's function, where N(x, x',t) samples the number of electrons, or equivalently, the number of noncrossing 
world lines in the interval x — x'. In the spin-incoherent regime, here realized by first taking J — > and then T — ► 0, 
the spin configurations all carry the same weight, and the probability that the world lines in the interval x — x' 
all have the same spin (as required for a non-zero contribution to the low-energy Green's function) becomes equal 
to 2~ N ( X ' X ' T \ Given this, the calculation of the Green's function G(x,x',t) translates into the calculation of four 
correlators for spinless time-boosted fermions, with the spin averaging factor 2~ N ^ X,X ,r ) properly inserted^, 

G(x,x',t) = (2- n ^ x - x '' t ^(x 1 t)¥(x',0)) 

= (2- n ( x ' x '> t ^„{x,t)^_{x' 7 0)) + (2- n( - x > x '' t ^ + {x,t)iPI(x',0)) 

+ {2- n ^ x - x '^^(x,t)^ + (x',0)) + (2- n ^ x - x '^^ + (x,t)^_(x',0)). (16) 

The operator N(x, x',t) can be expressed as 

N(x, x', t) = n\x - x'\ + 4=fo(*, t) - <p(x',0)}, (17) 

where n is the average electron (or world line) density, and with (l/y^)[ip(x, r) — ip(x',0)] — [cosh(i9) + 
sinh(^)] [<t>- (x, r) — </>_(— x,t) — 0_(a;',O) + x', 0)]/v / 7r the fluctuation term. 

The correlators in Eq. (TTBJ) are straightforwardly calculated by first using Eq. (fT3"]) to replace all occurrences of 
right-moving fermion fields by left movers, and then applying the bosonization formula (|14|) . Introducing relative and 
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center-of-mass coordinates, r = x — x' and R = (x + x')/2, respectively, we obtain for the G piece of the Green's 

function: 

G — (x, a;', r) = {2~ N ^ x '^_{x, r)^_ {x\ 0)) 

1 _ e -Qn2/n)~k F \r\ e ~i~k F r e iC- 1 1 



where 



2ira (a signr + v c t + ir) ((a signr + v c t) 2 + r 2 ) 2A1+2A2 

/ , \ 2A1+2A3 

/ ^{a 2 + {2R + r) 2 ){a 2 + {2R-r) 2 ) \ 

\ (a signr + v c t) 2 + AR 2 J ' 



ln2 / (a 2 + (2i? + r) 2 )(a 2 + (2i?-r) 2 ) ol (a signr + v c t) 2 

C — = -; — A In - : -5 h21n ^ 

47r V ((« sl S nr + v ct) 2 + AR 2 ) 2 a 2 



(18) 



ln2 / (a-i(r + 2.R))(Q;-i(r-2 J R)) ^ ^ ^ a signr + i; e r + ir \ ( 
47r V (a + z(r + 2i?))(a + i(r-2i?)) 11 a signr + w c r - ir / ' 

and with the exponents given by A x = f (^) 2 , A 2 = |(y + if - 2) , A 3 = \ (± - K) 
The G 1- part of the Green's function reads 

G-+(x,x',r) ee (2- w ^V_(x,r)V^',0)> 

1 _ e -(ln2/7T)fc I? |r| e -i2fc f .fl e iC_ + . 1 



2na (a signr + w c r + 2iR) ((a signr + v c t) 2 + 4R2)2Ai+2A 2 

\ 2Ai 

m signr + w c r) + r ) J 

2A 3 

: o o I . (20) 

(a signr + v c T) + r z I 



V(a 2 + (2i? + r) 2 )(a 2 4 


- (2R- 


r) 2 )(0 


y/(a 2 + (2R + r) 2 )(a 2 - 


f (2ii- 


. r )2)\ 



where 



la2^ a 2 + (2i? + r) 2 
-iv m ■ 



b + 4tt a 2 + (2R - r) 2 

In 2 / (a-«(r + 2iZ))(a-i(r-2.R)) a signr + u c r + ir\ 

+ -In — ; 77 + 2 m . 21 

4vr V (a + i(r + 2i?))(a + «(r - 2R)) a signr + v c t - ir J 

The expressions for the G ++ and G H pieces of the Green's function are immediately obtained from Eqs. (fT8|) - (f2"Tj) 

by using that 

C++ re', r) = G*__(x,x', r), 

G+^a/.r) - Gl + (ar ) ar',r). (22) 

Having derived the full expression for the boundary Green's function we can study boundary and bulk regimes by 
taking the proper limits. The bulk regime is defined by R 3> w c r, r. For the chiral G piece we obtain 



G (x, x' , r) = (2- n ( x > x ' t ^_(x,t)^_(x' ) 0)) 

1 

2na 



' -(ln2/7r)fcjr|r| — ifcpr iC- 



1 1 

x (^23) 

(a signr + v c t + ir) ((a signr + v c t) 2 + r 2 )~ 2Al+2A2 ' 

in agreement with the result for the bulk Green's function derived in Refs. [|| and ||. In contrast to the chiral part 
G (G++) of the Green's function, G \- (G-| ) decays with a power law and vanishes in the bulk limit, as it must: 

G- + (x ) a> ) T)~-^. (24) 

For the boundary case, defined by v c t >> r,R, the chiral and nonchiral parts of the Green's function show the 
same behavior, with the asymptotic scaling 

G — (x,x',t) ~ G-+{x,x',t) yj^- ( 25 ) 

This agrees with the result in Ref. [T3| (see also Ref. [Hj]). Importantly, our general result allows for a study of the 
crossover between boundary and bulk regimes. 
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IV. STM RESPONSE IN THE SPIN-INCOHERENT REGIME 



In this section we discuss a possible experimental probe of the spin-incoherent regime, using scanning tunneling 
microscopy (STM). Assuming that the STM tip couples only to the conduction electrons, the differential tunneling 
conductance dI(V,x)/dV is directly proportional to the local electron tunneling density of states (LDOS) p(V,x), 
where V is the applied voltage, and x is the position of the tip as measured from the boundary -2. The boundary may 
here be one of the endpoints of a quantum wire or a metallic nanotube, or generated dynamically^ by an impurity 
in the system (an antidot, or a deep "core level" hole created by an X-ray photon^). The LDOS is related to the 
boundary Green's function in Eq. (| 16(1 by 

p{x,un) fa -—lmG R {x,u)), (26) 

where the retarded Green's function Gr{x,uj) is obtained from the Fourier transform G(x,u> n ) of G(x, x, r) in Eq. 
(j 1 6(1 by analytically continuing iuj n — > lu + The ~ sign in Eq. (|26|) is a reminder that the expression for the 
Green's function in Eq. (116|) as derived above becomes exact only in the long-wavelength limit [cf. the text after 
Eq. ([5])]. The LDOS has a uniform [rapidly oscillating] part originating from the chiral [nonchiral Friedel-type] terms 
G — /++ (x,x,t) [G_ +/+ _(x,x,t)] in G(x,x,t), 

p(x, u) = p mli (x, lu) + cos(2k F x)p osc (x, lu). (27) 




FIG. 1: Uniform part of the LDOS, p un i, as a function of energy to for different boundary-to-tip distances x. Oscillations for 
intermediate distances are a result of interference of incoming and reflected charge modes at the boundary. 

Figure 1 shows p un { for three choices of distance from the boundary, obtained from Eq. (|16[) by numerically 
computing the integrals which define the analytically continued Fourier transforms. The most interesting feature is 
the oscillation pattern seen for "intermediate" distances where lux ~ v c (i.e., the crossover regime between boundary 
and bulk behavior). This pattern is due to the interference of incoming and reflected charge modes at the boundary. 
The oscillations, of wavelength 7tv c /x, give immediate information about the speed v c of the charge modes, and 
therefore about the strength of the electron-electron interaction. This unique signal of the spin-incoherent phase 
should be readily obtained in STM experiments as an oscillation pattern in the differential tunneling conductance as 
a function of the applied voltage V (which fixes the frequency for which the LDOS is probed). Let us here point out 
that the oscillations are present at any distance from the boundary, but that in the boundary [bulk] regime, defined 
by lux <C v c [lux 3> v c ] , the wavelength may become too long [short] to be possible to resolve experimentally. (Cf. Fig. 
1 where the wavelength for the oscillation for the case kpx — n is roughly twice larger than the displayed frequency 
interval.) It is important to stress that the Friedel terms [the last two terms in Eq. p^|) ] do not change the result 
qualitatively. This is in contrast to an ordinary Luttinger liquid where the interference of charge and spin modes 
with different speeds produce an oscillation pattern with very different properties^. In the case of a spin-incoherent 
system there is only one propagating charge mode with a single speed v c , and the oscillation pattern of p sc will be 
the same as that of p un i- 

Let us conclude this section by mentioning that other experimental probes of spin-incoherence have been suggested 
in the literature, including momentum resolved tunneling experiments on cleaved-edge overgrowth quantum wire s 15 ' 16 , 
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and quantum interference experiments on specially designed devices^. It would be very interesting to explore how the 
characteristic crossover behavior of the single-electron Green's function identified above may influence these proposed 
experiments. 



V. CONCLUSIONS 



We have calculated the exact single-electron boundary Green's function for one-dimensional spin-incoherent elec- 
trons in the low-energy, long- wavelength limit, using a bosonization scheme valid in the strong-coupling regime. The 
Green's function thus obtained correctly reproduces known results in the bulk and extreme boundary regimes. As 
revealed by the expressions in Eqs. (fT8| and (|20|) . the charge sector power-law scaling of the Green's function is 
highly sensitive to the presence of the boundary also at intermediate distances away from it, where the center-of-mass 
coordinate R ~ 0(v c t). This feature will strongly influence the local tunneling density of states, and may be probed 
experimentally by scanning tunneling microscopy of one-dimensional conductors of spin-incoherent electrons. 
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